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ABSTRACT 

The violent hierarchical nature of the A-Cold Dark Matter cosmology poses serious difficulties for the 
formation of disk galaxies. To help resolve these issues, we describe a new, merger-driven scenario 
for the cosmological formation of disk galaxies at high redshifts that supplements the standard model 
based on dissipational collapse. In this picture, large gaseous disks may be produced from high- 
angular momentum mergers of systems that are gas-dominated, i.e. M^as/{M.gs.s + Ad*) > 0.5 at 
the height of the merger. Pressurization from the multiphase structure of the interstellar medium 
prevents the complete conversion of gas into stars during the merger, and if enough gas remains 
to form a disk, the remnant eventually resembles a disk galaxy. We perform numerical simulations 
of galaxy mergers to study how supernovae feedback strength, supermassive black hole growth and 
feedback, progenitor gas fraction, merger mass-ratio, and orbital geometry impact the formation of 
remnant disks. We find that disks can build angular momentum through mergers and the degree of 
rotational support of the baryons in the merger remnant is primarily related to feedback processes 
associated with star formation. Nearly every simulated gas-rich merger remnant contains rapidly- 
rotating stellar substructure, while disk-dominated remnants are restricted to form in mergers that 
are gas-dominated at the time of final coalescence. Typically, gas-dominated mergers require extreme 
progenitor gas fractions (/gas > 0.8). We also show that the formation of rotationally-supported stellar 
systems in mergers is not restricted to idealized orbits, and both gas-rich major and minor mergers 
can produce disk-dominated stellar remnants. We suggest that the hierarchical nature of the A-Cold 
Dark Matter cosmology and the physics of the interstellar gas may act together to form spiral galaxies 
by building the angular momentum of disks through early, gas-dominated mergers. Our proposed 
scenario may be especially important for galaxy formation at high redshifts, where gas-dominated 
mergers are believed to be more common than in the local Universe. 

Subject headings: galaxies: formation - galaxies: evolution 


1. INTRODUCTION 

The conventional theory for the origin of disk galaxies 
in Cold Dark Matter (CDM) cosmologies involves the dis¬ 
sipational collapse of gas inside relaxed da rk matter halos 
formed through hierarchical clustering l|White fc Reefl 
11 97Stlffliimentlial et a,].l1984li . This theory endows form¬ 
ing disk galaxies with kinem atic properties inherite d 
from the halos that host them llFall fc Efstathioiilll98nll . 
Early simulations of cosm ological structure formation 
llBarnes fc EfstathioiJ Il9^ verified that during linear 
growth, the angular momen tum of halos g rows accord¬ 
ing to perturbation theory l|Peebleslll969ji . and by in¬ 
corporating these results subsequent work explained the 
flat rotation curves and sizes of disk galaxies through 
the combination of gas dissipation, halo spin, and the 
adiabatic response of dark matter to baryonic infall 
llBlumenthal et al.ll98(Tl:|Mo et al.ll998fl . The initial suc- 
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cess of the dissipational collapse model led to its estab¬ 
lishment as the favored theory for the formation of disk 
galaxies. However, numerical modeling revealed that the 
theory was incomplete. 

While smooth, dissipational collapse may be rele¬ 
vant to disk galaxy formation in quiet environments, 
simulations and semi-analytic modeling have demon¬ 
strated that it is rare for halos to acquire most 
of their angular momentum via the quiescent accre¬ 
tion of tidally torqued material M itvitska et af] 120021 
iMaller fc Deke]lE002HD’Onghia fc BuAe^2004h . Simu¬ 

lation and observational work has shown that the spe¬ 
cific angular momentum content of dark matter ha¬ 
los differs from exponential dis ks (|Bu1]ock et aUl^flflll 
Ivan den Bosch et alJl200lt 1200211 as halos contain com¬ 
paratively more low angular momentum material that 
might produce bulges or overly c entrally concentrated 
disks llvan den Boschl IT998i 1200111 . Furthermore, sim¬ 
ulations of mass accretion onto galaxies indicate that 
the hierarchical nature of the ACDM cosmology leads 
to the destr uction of disks when dissipative effects 
are neglected (iTo th fc Ostrikeilll992l: iQuin n et aUlToO.Jl: 

iWalker et al.lll990HVelazauez fc Whitelll999ll. Models of 

galaxy collisions including dissi pation (iHernaui^ 119891: 

iBarnes fc Hernauistlll 99 111 199011 show that gas can loose 
angular momentum owing to gravitational torques dur¬ 
ing these events. In cases when the interstellar medium 
(ISM) is isothermal and relatively cold and the gas frac- 
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tion of the galaxies is small (~ 10% of the baryons), the 
subsequent inflow of gas into the centers of the merger 
remnants forms a roughly sp herical stellar distribution 
throu gh a luminous starburst llMihos fc Hernauistll994cL 
[THi), leaving objects that have essentially no extended 
stellar disks. Cosmological simulations of disk formation 
with similar physics produce galaxies that are too cen¬ 
trally concentrated, form too many stars, contain overly 
dominant bulges, and l ack angular momentum compared 
with nearby spirals llKatz fc CunnI ITo^ iKatz et al.l 
1992t iNavarro fc Whitelll994 iSteinmetz fc Mulleilll995l: 

Navarro fc Steinmetzll2000(). 

Efforts to resolve these discrepancies in a cosmolog¬ 
ical context have mainly relied on in creasing the im¬ 
pact of feedback from star forma tion tWeiLet_aI 
Somnier;_Larsene^^yil999L|2003jThacke^^_Coudmi^ 

200^]20^||^Aadr^t^^]2003a ^Cpye^^^^^^ liiol 

Robertson et al 12004 lOkamoto d al I2005D . Cosmolog¬ 
ical simulations that utilize strong feedback from star 
formation in the form of a pressurized ISM to produce 
disk galaxies either produce too few exponential d isks 
compared with observations llRobertson et alj|2004l or 
yield ga laxies that still con tain a significant bulge com¬ 
ponent llAbadi et alJl2003^ . While these recent simula¬ 
tions have had some success, taken together the results 
suggest that the theory of structure formation remains 
faced with an apparent incompatibility between models 
of disk formation and the hierarchical growth of struc¬ 
ture. 

To help resolve this incompatibility, we propose a new, 
merger-driven theory for the formation of disk galax¬ 
ies at high redshifts that supplements standard dissipa- 
tional collapse. In our scenario, gas-dominated merg¬ 
ers, where the interacting systems have a gas frac¬ 
tion /gas = Mgas/(Mgas + Al*) > 0.5 at the height 
of the merger, allow for the formation of rotationally- 
supported disks in remnants if energetic feedback mech¬ 
anisms limit the conversion of gas into stars. Gas pressur¬ 
ization f rom a multiphase ISM by feedback from star for¬ 
matio n ijMcKee fc Ostrik^ 119771: iSnringel fc HernauistI 
l2?inl stabilizes gas disks and allows mergers between 
gas-dominated galaxies to produce large, smoothly dis¬ 
tributed, rapidly rotating stellar disks. Our proposed 
scenario builds upon the results of previ ous simu l ations 
of r emnant disk formation in me rgers byE arne i il2nnl 
and lSoringel & HernauistI ll200,5tl . 

The merger-driven scenario for disk formation transi¬ 
tions naturally into the mer ger-driven sceii ario for ellip¬ 
tical galaxy formation (e.g. iToomr el 1197111 as star for¬ 
mation in galaxies over time reduces the gas content of 
progenitors and prevents the further occurrence of merg¬ 
ers involving gas-dominated systems. Mergers between 
gas-rich disk galaxies at intermediate redshifts (z ft; 1 — 3) 
are still expected to produc e realistic elliptical ga lax- 
ies nearly universally (see e.g. lRobertson et 'ml200(i^ . as 
prior quiescent star formation and interactions will de¬ 
crease their gas content below the necessary /gas ^0.5 
needed to produce a disk-dominated remnant. 

Below we describe the simulation methodology used to 
demonstrate the merger-driven formation of disk galax¬ 
ies (©, incl uding our methods for st udyi ng the effects 
o f fee dback (EH, orbital geometry (E3i, gas fraction 
(E3|), merger mass-ratio (El, and progenitor mass 
(El on remnant properties. Our method for analyzing 


the structural and kinematic properties of the remnants 
is detailed in il2.fil a,nd the results of the numerical exper¬ 
iments are presented in E We discuss and summarize 
our results in E and E 

2. METHODOLOGY 

To demonstrate the feasibility of our scenario for 
merger-driven disk formation, we have performed a suite 
of merger simulations of isolated, gas-rich disk galax¬ 
ies, with and without feedback from accretion onto su- 
permassive black holes (BHs). The sim ulations use 
Smoothed Particle Hydrody namics (SPH) l|Lucvl IT977t 
iGingold fc MonaghanI I 197 D to evol ve the gas and in¬ 
corporate a multiphase ISM model tV&Des_et_al 1997; 
IRjii^nignjj_j^aTg^g 3 ^ iSoringel fc Hernauist 200,1 ; 
Snringel et abl l2005hj) where we control the pressur¬ 
ization (/EOS of star-forming gas by interpolating be¬ 
tween isothermal gas (gsos = 0) with a temperature 
Teg = 10"* K and a stiffer multiphase medium (gsos = 
1) with a higher effective temperature Teg « 10^ K 
by varying the equa tion of state (EOS; for details see 
ISoringel et al.ll^05bH . Star formation in the dense ISM 
(E > IOM 0 p c~^~) occurs on a timescale chosen to match 
observations Jk ennicu 13 mn). BH growth is modeled 
ISoringel et al.ll2005bH by spherical accretion onto BH 
“sink” particles, with a fraction Cf of the released grav¬ 
itational potential energy of the accreted gas being de¬ 
posited as thermal feedback into the gas surrounding the 
BH. For the simulations and accretion model considered 
here, Cf = TythermCr = 0.005, where Cr « 0.1 is the frac¬ 
tion of rest mass energy radiated during the growth of the 
BH (fixed by standard accretion theory) and the thermal 
coupling efficiency J/therm = 0.05 is chosen ter repr oduce 
the local Mbh — o' relation llDi Matteo et al.ir2005D . 

We consi der two merger progenitors, a Vvir = 160 km 
system llSoringel fc Hernauistll200,^ with a rotation 
curve and virial mass = 9.5 x IO^^/i^^Mq similar to 
the Milky Way and a smaller, Vvir = 80 km s“^ system 
with virial mass M^g = 1.2 x The galaxies 

are comprised of exponential gas and stellar diste embed¬ 
ded in dark matter halos with a lHernauistl (1199011 profile 
having scale lengths corresp onding to a Navarro-Frenk- 
White l|Navarro et al.lll9^ concentration of cnfw = 9. 
The progenitors each contain Ng^s = 40000 gas, IVgtars = 
40000 stellar, and Ahm = 180000 dark matter halo parti¬ 
cles, and may also include a BH seed of lQ^h~^MQ that is 
allowed to accrete in the manner described above. Once 
the BH particles fall within our force resolution limit, the 
sink particles combine into a single supermassive BH. 

To characterize the conditions under which our pro¬ 
posed scenario for merger-driven disk formation is most 
successful, we simulate a variety of merger scenarios vary¬ 
ing the ISM pressurization ((/eos), gas fraction (/gas), 
merger mass-ratio, disk orientations, and orbits. The 
complete suite of simulations is described below, and a 
summary is provided in Table ^ For convenience we 
will refer to each model with a two letter designation, 
with the first letter signifying the gas fraction /gas and 
the second letter noting the ISM pressurization gEOS- 
When necessary, we add suffixes to describe e.g. orbital 
or merger mass-ratio variations. The simulation desig¬ 
nations are detailed in the following sections. 

2 .1. Feedback Study 
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Table 1 . Merger Models 


Model 

^ ir , A 
[km s“^] 

[km S“^] 

/gas 

IJEOS 

^peri 

[h~^ kpc] 

Hnit 

[h~^ kpc] 

9 i 

[degrees] 

4^1 

[degrees] 

02 

[degrees] 

4>2 

[degrees] 

BH 

GA 

160 

160 

0.99 

0.1 

6.0 

70.0 

0 

0 

0 

0 

Y 

GB 

160 

160 

0.99 

0.25 

6.0 

70.0 

0 

0 

0 

0 

Y 

GG 

160 

160 

0.99 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

GD 

160 

160 

0.99 

0.75 

6.0 

70.0 

0 

0 

0 

0 

Y 

GE 

160 

160 

0.99 

1.0 

6.0 

70.0 

0 

0 

0 

0 

Y 

No BHs 

GAn 

160 

160 

0.99 

0.1 

6.0 

70.0 

0 

0 

0 

0 

N 

GBn 

160 

160 

0.99 

0.25 

6.0 

70.0 

0 

0 

0 

0 

N 

GGn 

160 

160 

0.99 

0.5 

6.0 

70.0 

0 

0 

0 

0 

N 

GDn 

160 

160 

0.99 

0.75 

6.0 

70.0 

0 

0 

0 

0 

N 

GEn 

160 

160 

0.99 

1.0 

6.0 

70.0 

0 

0 

0 

0 

N 






^peri 

Variations 






GCrA 

160 

160 

0.99 

0.5 

2.7 

70.0 

0 

0 

0 

0 

Y 

GGrB 

160 

160 

0.99 

0.5 

16.0 

70.0 

0 

0 

0 

0 

Y 

GCrC 

160 

160 

0.99 

0.5 

32.0 

70.0 

0 

0 

0 

0 

Y 

GGrD 

160 

160 

0.99 

0.5 

48.0 

140.0 

0 

0 

0 

0 

Y 

Orbital Variations 

GGoC 

160 

160 

0.99 

0.5 

6.0 

70.0 

180 

0 

180 

0 

Y 

GGoE 

160 

160 

0.99 

0.5 

6.0 

70.0 

30 

60 

-30 

45 

Y 

GCoF 

160 

160 

0.99 

0.5 

6.0 

70.0 

60 

60 

150 

0 

Y 

GGoO 

160 

160 

0.99 

0.5 

6.0 

70.0 

-109 

30 

71 

-30 

Y 

GGoP 

160 

160 

0.99 

0.5 

6.0 

70.0 

-109 

90 

180 

0 

Y 






/gas 

Variations 






DG 

160 

160 

0.4 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

EC 

160 

160 

0.6 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

EC 

160 

160 

0.8 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

Minor Mergers 

DCm 

160 

80 

0.4 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

ECm 

160 

80 

0.6 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

FCm 

160 

80 

0.8 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

= 80 km s ^ Progenitors 

DCs 

80 

80 

0.4 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

ECs 

80 

80 

0.6 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 

FCs 

80 

80 

0.8 

0.5 

6.0 

70.0 

0 

0 

0 

0 

Y 


To study the impact of energetic feedback mechanisms 
on the formation of disks in gas-rich mergers, we per¬ 
form 12 simulations of major mergers between Vvir = 160 
km s“^ galaxies with disks consisting almost entirely of 
gas (/gas = 0.99, “G” runs). The ISM pressurization 
is incrementally increased over the values Qeos = O-lj 
0.25, 0.5, 0.75, and 1.0, stiffening the equation of state 
of the ISM and increasing the sound speed of the gas at 
a given density and temperature. The simulations are 
referred to as models GA-GE in Table ^ We perform 
each simulation twice, once with and once without the 
inclusion of growing supermassive BHs in the progenitor 
systems. The models without BHs are listed with the 
suffix “n” in Tabled and the model GGn (/gas = 0.99, 
qF.ns = 0.5) corresponds to the simulation performed by 
iSpringel fc Hernouisd lj2005jl . 

The galaxies merge on coplanar, prograde orbits with 
an initial separation nnit = 70h~^ kpc and pericentric 
distance Cperi = 6h~^ kpc. The orbit for the feedback 
study r uns was chosen to reproduce the encounter simu¬ 
lated bv lSoringel fc HernauistI ll2005ll . By using the same 


orbit for these simulations, we explore the importance of 
feedback for the formation of disks in mergers separate 
from the role of angular momentum. We explore the 
role of angular momentum in producing remnant disks 
in mergers through additional simulations described in 

a 

2.2. Orbital Study 

The original orbital and disk angular momentum of the 
progenitors will affect the ability of disks to form from 
mergers as angular momentum present in residual gas left 
over from the merger can lead to rapidly rotating struc¬ 
tures in the remnant. The formation of remnant stellar 
disks have been demo nstrated in interactions with pro¬ 
grade, coplanar orbits l|Springel fc Hernauistll200,^ . and 
in simulations without star formation gaseous remnant 
disks have been shown to form in me rgers with both po¬ 
lar and inclined orbits lIBarne iEool. 

To judge the importance of orbital and disk angular 
momentum on a merger-driven scenario for disk galaxy 
formation, we perform a study of 9 additional Simula- 
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tions. In four models, we vary the pericentric passage 
distance rperi of the merger to alter the orbital angular 
momentum of the interaction. Typically, increasing rperi 
will increase the angular momentum of the orbit and lead 
to larger, more rotationally-supported disks. However, 
mergers with larger rperi undergo inefhcient dynamical 
friction, avoid strong angular momentum loss during the 
first passage, and therefore take longer to merge. Quies¬ 
cent star formation will reduce the effective gas fraction 
of the progenitors before the disks merge in these cases, 
and may lead to larger stellar spheroids. These two com¬ 
peting effects will be present in cosmological environs to 
some degree, and for this reason we do not attempt to 
adjust for the decrease in gas fraction in wide-orbit merg¬ 
ers. In addition to a coplanar orbit with Cperi = 
kpc, which is approximately twice the progenitor disk 
scale length R^, we simulate one model with a smaller 
pericentric passage distance (rperi = Rd, model GCrA 
in Table P) and three with larger pericentric passage 
distances (rperi = 0.1i?vir, 0.2i?vir, and 0.4i?vir, models 
GCrB-GCrD in Tabled]). 

In the remaining hve models of the orbital study, we 
vary the orientation of the progenitor disks. While vary¬ 
ing the pericentric passage distance increases the orbital 
angular momentum, much of this angular momentum 
must be lost to the dark matter halo through dynam¬ 
ical friction before the disks can merge. The original 
disk angular momenta of the progenitors may then sub¬ 
stantially add to the angular momentum of the rem- 
na nt disk. In addition to the c oplanar orbit considered 
bv lSoringel fc HernauistI (|200^ w e selec t five additional 
disk orientations from the^^me mui) study of equal- 
mass mergers of stellar disk galaxies. These orientations 
range from retrograde-retrograde coplanar {9i = 180, 
4’! = 0, 02 = 180, (t >2 = 0, model GCoC) to prograde- 
retrograde polar (01 = 60, (j)i = 60, 02 = 150, 4>2 = 0, 
model GGoF), with a complete list provided in Tabled 
The results of these simulations are presented in an 

2.3. Gas Fraction Study 

The production of remnant disks in mergers will de¬ 
pend on the gas fraction /gas of the progenitor disk sys¬ 
tems. Stellar disks present before the merger will be 
substantially heated during the merger and form stellar 
spheroids. Only galaxies with signihcant gas fractions 
before the merger will be able to re-form gaseous disks 
after the galaxies collide. To characterize the progenitor 
gas fraction necessary to form a substantial remnant disk, 
we perform three simulations with growing BHs, ISM 
pressurization (Jeos = 0.5 and gas fractions /gas = 0.4, 
0.6, and 0.8 (models DC, EC, and FC in Tabled). The 
results of these simulations are presented in a 

2.4. Minor Merger Study 

In CDM cosmologies, the frequency of galaxy mergers 
typically inc reases with the mas s ratio of the merging 
pair (see e.g. lLacev fc Golelll99.3ll . As mentioned above, 
previous theoretical studies have suggested that dissipa¬ 
tionless interactions between satellites and disk galaxies 
can destroy disks even if the satellite mass is relatively 
small (~ a few percent of the disk mass, e.g. lQuinn et abl 
rr^ . Minor mergers with progenitor mass ratios of 
3 : 1 have been shown to produce s mall gaseous d isks 
in simulations without star formation llBarne i|2002). In 


simulations including gas dissipation and star formation, 
min or mergers have been show n to induce starbursts 
fe.g iMihos fc HernauistI ITQQdbll . Cosmological simula¬ 
tions suggest that some minor mergers may also add 
mass to the thick component of disk galaxy systems 
ijAhadi et a,1.ll^fl3hfl . Minor mergers may then greatly 
influence the formation and survival of disk galaxies 
even though the mass accretion rate from minor merg¬ 
ers is typically smaller than that from major mergers. 
To address the relevance of minor mergers for remnant 
disks, we rerun the gas fraction study simulations from 
El substituting a I4ir = 80 km s ^ galaxy for the 
second progenitor. These mergers with mass ratios of 
TO 2 /TO 1 ~ 8 : 1 will help characterize the influence of 
minor mergers on remnant disk properties and provide 
examples of models that approximate cosmologically fre¬ 
quent merger events. The results of these simulations are 
presented in El 

2.5. Small Progenitor Study 

The scale dependent physics of gas dissipation and star 
formation may influence the ability of remnant disks 
to form in mergers. Gas cooling is typically more effi¬ 
cient in smaller mass halos as the virial temperature of 
those systems approach the regime where hydrogen and 
helium recombination and collisional io nization enables 
large cooling rates (see e.g. H Gas densities 

in the disks of low mass systems are typically lower than 
in more massive systems, leading to longer gas consump¬ 
tion timescales. Quiescent star formation in low mass 
systems during the phases of the merger when the sep¬ 
aration between disks is large and gravitational torques 
are weak will be lower than for high mass systems, which 
may correspond to a larger gas fraction in the progenitor 
disks when the systems actually merge llRobertsoTi et alJ 

Hoi). 

To gauge the net impact of these effects on the forma¬ 
tion of remnant disks in mergers of different mass scales, 
we perform three additional simulations of equal mass 
mergers between T4ir = 80 km s“^ galaxies wit gas frac¬ 
tions of /gas = 0.4, 0.6, and 0.8 for the progenitor sys¬ 
tems (models DCs-FGs in Ta ble P ). The results of these 
simulations are presented in El 

2 .6. Analysis 

Each simulation is evolved over the timescale needed 
for the progenitor systems to merge, corresponding to 
T = 3.0 — 10.3 Gyr depending on the orbit. To illustrate 
the structure of the merging system during the collision, 
Figure n shows the time evolution of the gas in merg¬ 
ers with a strongly pressurized (^eos = 1-0) ISM with 
(model GE, bottom panel) and without (model GEn, 
top panel) feedback from accretion onto supermassive 
BHs. Each panel shows a 140 kpc x 140 kpc square area 
around the center of mass of the interacting system. At 
early times in the merger (T = 0.35 Gyr, far-left panels), 
tidal interactions between the galaxies distort them, but 
as signihcant accretion onto the BHs has yet to occur, the 
overall gas distribution remains similar between the sim¬ 
ulations. As the gas in the simulation with BHs (bottom 
panels) begins to experience thermal feedback from BH 
accretion, a hot halo of diffuse, low-angular momentum 
gas expelled through a wind surrounds the interacting 
galaxies. By the end of the simulations at T = 2.0 Gyr 
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T-a.SSGyr T-OTIGyr T-1.4 £Cvt T-2.DDCvr 

GE 



T-a.SSGyr T-OTI^yr T-1.4 £Gvt T-2.DDGvr 


Fig. 1.— Merger of two gas-rich disk galaxies, with (bottom, model GE) and without (top, model GEn) supermassive black hole accretion 
and feedback. Shown is the gas surface density of the galaxies in each 140 X 140 kpc square panel, demonstrating that the high-angular 
momentum merger results in a remnant disk. The effects of black hole feedback can be seen over the 2 Gyr timescale of the merger as a 
diffuse, hot halo and the lowered central gas density in the remnant (bottom panel). 


(far-right panels), rotating gas disks remain, eventually 
forming stellar disks. The remnant containing a BH has 
cleared out some gas from the inner regions of the disk, 
reducing the central stellar density and producing a dif¬ 
fuse hot gas envelope. 

After the merger has completed, the structural and 
kinematic properties of the remnant are measured from 
the simulation data. The structural components of the 
merger remnants, which might include a bulge, a thin 
or thick disk, or an extended spheroid, are comprised of 
stars that form at characteristically different times dur¬ 
ing the interaction. The dynamically hot components of 
the remnant typically form before the final coalescence of 
the merger, while the stars in the colder components of 
the remnant mostly form after the height of the merger 
when the role of tidal forces has decreased. The star 
formation rate (SFR) of the merging system (see Figure 
13 upper left panel) can be used define selection crite¬ 
ria for stars likely to populate separate dynamical com¬ 
ponents of the remnants based on their formation age. 
The peak in the SFR that occurs during the height of 
the merger can be approximated by an exponential rise 
and decay. In what follows, we identify “old” stars in 
the remnants as stellar particles forming earlier than be¬ 
fore three e-folding times before the center of the SFR 
peak and “new” stars as those forming later than three 
e-folding times after the peak. 

The old and new stars in the remnant are then used to 
measure the separate kinematics of the dynamically hot 
and cold stellar structures of the remnant. Typically, the 
old stars in the remnant have relatively low rotational ve¬ 
locities (Figure 13 lower left panel) and are not rotation- 
ally supported with the ratio of rotational velocity to ve¬ 
locity dispersion I^ot/o’ < 1 (Figure[3 lower right panel). 
New stars formed in the remnant after the height of the 
merger typically have high rotational velocities and are 
rotationally supported. Stars formed during the burst 


can have kinematic properties intermediate between the 
old and new stars, and therefore are not analyzed as a 
separate stellar component. 

The mass of the new stellar component forming in the 
remnant disk can be comparable or large than the old, 
dynamically heated bulge or extended spheroid. The 
mass of the stellar structures in the remnant can be 
measured from the surface mass density profile. After 
the merger the stellar angular momentum of the rem¬ 
nant is measured and used to define a possible “disk”- 
plane perpendicular to the rotation axis of the remnant. 
The plane defined normal to the stellar angular momen¬ 
tum vector captures the mean stellar rotation in the rem¬ 
nants we study, and only underestimates rotation when 
rotationally-supported stellar structures are misaligned 
with the total stellar angular momentum. The stellar 
surface mass density of each remnant is then measured 
in logarithmically-spaced annuli and is modeled by a 
multi-component system (Figure [3 upper right panel) 
as a function of radius r that may include an exponen¬ 
tial bulge 


Sb(r) = 

(1) 

an exponential thin disk 


Ed(r) = 

(2) 

an exponential thick disk 


Ed,thick(r) = 

(3) 


and a n extended stellar spheroid with a ide Vaiicouleurj 
l|1948D profile 


Es(r) = (4) 

Each remnant is fitted with individual thin disk, bulge- 
thin disk, bulge-thin disk-thick disk, bulge-thin disk- 
spheroid, and bulge-thin disk-thick disk-spheroid mod¬ 
els with the requirements that the bulge scale length 
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Fig. 2. — Example analysis performed for each merger simulation, shown for an equal mass merger between Vvir = 160 km s“^ galaxies 
with gas-rich disks (/gas = 0.99), strongly pressurized ISM (^eos = l-O), and growing supermassive black holes (model GE from Table \T\. 
The peak of the star formation rate (SFR) measured during the simulation is used to identify the height of the merger (upper left panel). 
An exponential rise and decay is fitted to the SFR peak (dotted line) and used to categorize the stellar content of each remnant based 
on its formation time. Stars forming before three e-folding times prior to the maximum of the SFR rate are labeled as “old” stars (left 
dashed line), while stars forming three e-folding times after the SFR peak are labeled as “new” stars (right dashed line). Typically the 
new stars forming from remnant disks have large rotational velocities (lower left panel) and are rotationally supported with the ratio of 
rotationally velocity to velocity dispersion Vrot/o" > 1 (lower right panel). The surface mass density S(r) (upper right panel) is modeled as 
a multi-component system consisting of an exponential bulge, think disk, thick disk, and extended spheroid. The model shown here has a 
surface mass density best modeled by a bulge-thin disk—thick disk system. 


ab < 1 kpc, the thick disk scale length i?d, thick is larger 
than the thin disk scale length i?d, and the spheroid scale 
length Os > Ob- A simple estimate of the success 
of each model in reproducing the surface mass density 
profile of the remnant is used to determine the relative 
mass in bulge, disk, and spheroid components. When the 
model fits suggest that either a thick disk or an extended 
spheroid models provide comparable fits to the projected 
stellar distribution, the rotational support Vot /cr oi the 
stars at large radii is used to select a preferred model. In 
this case, dispersion supported systems (K-ot/o" < 1) are 
assigned extended spheroids and rotationally-supported 
systems (Vot/c > 1) are assigned thick disk components. 
The best-fit surface mass density models are reported in 
Table 121 including the remnant bulge-disk-spheroid mass 
ratios. For comparison, Table|21also lists the best fit sur¬ 
face mass density parameters for a Kir = 160 km s“^ 
progenitor (model ICs) with a gas fraction of /gas = 0.4 
measured using the same method. 

3. RESULTS 


The suite of simulations provides an extensive set 
of remnants that characterize the ability of extremely 
gas-rich mergers to produce both disk components in 
spheroid-dominated systems and disk-dominated sys¬ 
tems whose structure and kinematics closely resemble 
large spiral galaxies. In nearly every gas-rich merger 
we simulated (save model DCs, see below) the new 
stars formed after the peak of the SFR produced a 
rotationally-supported structure. These remnant disks 
are ubiquitous in gas-rich mergers, but the extent to 
which they dominate the stellar content and the degree 
of rotational support shows trends with the strength of 
pressurization in the ISM, the presence of growing su¬ 
permassive BHs, progenitor gas fraction, orbit and disk 
orientation, and progenitor mass ratio. Below, we con¬ 
sider the impact of these merger properties in turn. 

3.1. Feedbaek Study 

The kinematic properties of the remnants exhibit sig¬ 
nificant quantitative differences as the pressure support 
of the ISM is varied. The stellar rotation curves of the 
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Fig. 3. — Stellar rotation curve of merger remnants, for models with differing ISM pressurization. Shown is the rotation of stars formed 
before the merger (open triangles), after the merger (open circles), and for all stars in the system (filled circles). The amount of rotation in 
the disk remnant correlates strongly with gas pressurization, increasing from the upper left panel (gsos = 0.1, model GAn) to the bottom 
middle panel ((Jeos = 1-0, model GEn). The nearly isothermal ISM model (qeos = 0-1, model GAn) produces almost no net rotation in 
the stellar disk, while moderately to strongly pressurized models (qeos = 0.5 — 1.0, models GCn-GEn) yield a rotating stellar component. 
The lower right panel shows a strongly pressurized ISM model (iJeos = 1-0, model GE) with BHs. 




r [h'' kpc] 


r [h'* kpc] 


Fig. 4.— Remnant stellar mass surface density with a strongly pressurized (gEOS = 1-5) ISM, with (right panel, model GE) and without 
(left panel, model GEn) black hole feedback. Shown is the measured stellar surface density (open circles), and bulge (dot-dashed line), thin 
disk (dashed line), thick disk (dotted_Jine2j_^dcomposite (solid line) mass models fits. We fitex£onential_^sks_to^the central bulge and 
stellar disk components and a|de_yaucouleu^ profile to the spheroid. As discussed by |Itobertson_et_^y (20061), black hole feedback 

reduces slightly reduces the mass of the central stellar spheroid and increases its effective radius. The relative sizes of the bulge and disk 
components are reported in Table 


merger remnant in the (?eos = 0.1,0.25,0.5, 0.75, and 1.0 
models without BHs are plotted in Figure 13 for a pro¬ 
genitor gas fraction of /gas = 0.99. The strength of the 
rotating stellar component correlates with the pressuriza¬ 
tion of the ISM. The least effective ISM model for pro¬ 
ducing a remnant disk is the weakly pressurized, nearly 
isothermal model ((?eos = 0.1, model GAn, upper left 
panel) that yields a slowly rotating, clumpy stellar rem¬ 


nant. The gaseo us progenitor d isks in model GAn are 
not stable against iToomrel 119641) instability as the sound 
speed of the weakly pressurized gas cannot support the 
massive gas disk against self-gravity. The gaseous disk 
efficiently converts into stellar clumps before and during 
the collision, producing a remnant with very little av¬ 
erage rotation (maxVrot ~ 90 km s“^). However, the 
scant amount of gas remaining after the final coalescence 
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Table 2. Merger Remnants 


Model 

Sb 


Sd 

Rd 

Ss 


^d,thick 

-^d,thick 

Mb 

Md 

Ms 

-^d,thick 

B:D:S 

ICs 



4.9e8 

2.73 






2.3el0 



0.0:1.0:0.0 

BHs 

GA 

6.9el0 

0.24 

2.1e8 

3.50 

8.3e7 

4.10 



2.6el0 

1.6el0 

3.1el0 


1.6:1.0:1.9 

GB 

3.9el0 

0.41 

2.1e8 

4.48 





4.2el0 

2.6el0 



1.5:1.0:0.0 

GC 

2.0ell 

0.17 

1.6e9 

1.55 



3.3e7 

5.19 

3.8el0 

2.5el0 


5.7e9 

1.2:1.0:0.0 

GD 

7.3el0 

0.18 

1.9e9 

1.42 

8.0e7 

2.89 



1.6el0 

2.4el0 

1.5el0 


1.0:1.6:1.0 

GE 

l.Sell 

0.18 

5.9e9 

0.65 



5.2e8 

2.26 

2.7el0 

1.6el0 


l.OelO 

1.0:1.2:0.0 

No BHs 

GAn 

4.3el0 

0.28 



1.6e8 

3.88 



2.1el0 


5.6el0 


1.0:0.0:2.5 

GBn 

5.0el0 

0.28 

8.2e7 

4.45 

3.7e8 

1.97 



2.6el0 

l.OelO 

3.3el0 


2.5:1.0:3.2 

GCn 

l.Sell 

0.18 

1.3e9 

1.61 

1.8e7 

4.15 



4.1el0 

2.1el0 

7.3e9 


5.6:2.8:1.0 

GDn 

1.9ell 

0.18 

1.2e9 

1.48 



1.4e8 

3.21 

4.1el0 

1.6el0 


9.6e9 

1.5:1.0:0.0 

GEn 

2.7ell 

0.12 

1.3el0 

0.45 



8.5e8 

2.02 

2.8el0 

1.7el0 


2.1el0 

1.0:1.4:0.0 

Orbital Study 

GCoC 

7.0el0 

0.25 

1.5e9 

1.78 



3.0e7 

5.39 

2.7el0 

3.0el0 


5.5e9 

1.0:1.3:0.0 

GCoE 

5.7el0 

0.27 

6.0e8 

2.14 

1.6e8 

2.30 



2.6el0 

1.7el0 

2.0el0 


1.5:1.0:1.1 

GCoF 

4.7el0 

0.25 

1.7e9 

1.96 





1.9el0 

4.2el0 



1.0:2.1:0.0 

GCoO 

7.2el0 

0.22 

1.5e9 

1.94 



9.3e6 

5.59 

2.3el0 

3.7el0 


1.8e9 

1.0:1.6:0.0 

GCoP 

6.2el0 

0.24 

1.6e9 

1.43 

9.0e7 

2.90 



2.3el0 

2.1el0 

1.7el0 


1.3:1.2:1.0 

GCrA 

5.7el0 

0.32 

2.0e9 

1.28 

2.8e7 

4.29 



3.7el0 

2.1el0 

l.lelO 


3.1:1.8:1.0 

GCrB 

2.4el0 

0.35 

8.6e8 

2.04 



6.1e7 

6.02 

1.9el0 

2.2el0 


1.3el0 

1.0:1.8:0.0 

GCrG 

4.0e9 

0.54 

1.5e8 

3.32 

l.OeS 

4.31 



7.5e9 

l.OelO 

4.4el0 


1.0:1.4:5.9 

GCrD 


0.67 

4.6e7 

5.57 

2.0e8 

3.42 




9.1e9 

5.4el0 


0.0:1.0:5.9 

Gas Fraction Study 

DC 

3.3el0 

0.25 

2.7e7 

6.15 

1.6e8 

3.65 



1.3el0 

6.6e9 

4.9el0 


2.0:1.0:7.5 

EC 

6.7el0 

0.21 

1.9e7 

4.71 

1.7e8 

3.38 



2.0el0 

2.7e9 

4.6el0 


7.6:1.0:17 

EC 

2.4ell 

0.11 

7.7e9 

0.54 

2.9e7 

4.89 

2.7e8 

2.77 

1.9el0 

1.4el0 

l.OelO 

1.3el0 

1.1:1.7:1.0 

Minor Merger Study 

DCm 

4.3e9 

0.57 

3.7e8 

3.36 





9.4e9 

2.6el0 



1.0:2.8:0.0 

ECm 

3.6e9 

0.58 

3.3e8 

3.14 

1.3e7 

4.70 



7.6e9 

2.0el0 

6.8e9 


1.1:3.0:1.0 

FCm 

7.5e9 

0.36 

3.7e8 

3.06 

1.3e7 

3.66 



6.3e9 

2.1el0 

4.2e9 


1.5:5.1:1.0 







V^ir 

= 80 km s 

^ Study 






DCs 

8.7e9 

0.12 

1.0e9 

0.33 

6.7e7 

2.15 



8.5e8 

7.2e8 

7.1e9 


1.1:1.0:9.8 

ECs 

3.3el0 

0.10 

9.1e8 

0.46 

3.5e7 

2.43 



2.1e9 

1.2e9 

4.8e9 


1.7:1.0:3.9 

FCs 

3.2el0 

0.10 

1.5e9 

0.34 

5.0e7 

1.92 



2.0e9 

1.2e9 

4.2e9 


1.7:1.0:3.4 


in this merger still manages to produce a very small, 
rapidly rotating disk. Increasing the ISM pressurization 
to gEOS = 0.25 (model GBn, upper middle panel) greatly 
reduces the presence of stellar clumps owing to instabili¬ 
ties and roughly doubles the radial extent of the rotating 
component, but does little to improve the average ro¬ 
tation of the remnant. Overall, low-pressurization ISM 
models do a poor job at producing remnants with sub¬ 
stantial average rotation but can produce small, rapidly 
rotating remnant disk components. 

Merger-driven disk formation will require a moderately 
to strongly pressurized ISM to provide enough kinematic 
stability to the gas to avoid forming too many stars be¬ 
fore the final merger. The moderately to strongly pres¬ 
surized models (gsos = 0-5 — 1-0) GCn-GEn, upper 
left-lower middle panels) produce remnants that contain 
rapidly rotating stellar disks in addition to bulge or ex¬ 
tended spheroid components. During the merger, the 
pressurization from the stiffer equation of state makes 
the ISM less compressible, reduces the central gas den¬ 
sity, and leads to less conversion of gas into stars. A 


larger fraction of stars in the final remnant form after 
the height of the merger, increasing the mass of the disk 
component and the average rotational velocity of the 
stars. Just before the final coalescence, the moderately 
and strongly pressurized models have /gas ~ 0.5 — 0.6, 
much larger than the /gas ~ O.I in the gEOS = 0.1 
model and /gas = 0.26 in the gEOS = 0.25 model. The 
abundance of gas present during the final coalescence 
in the moderately and strongly pressurized models al¬ 
lows for stellar disks to form from gaseous remnant disks 
(e.g. iSoringel fc Hernaiii^ l200,5ll . Figure 0 shows the 
stellar surface mass density of two highly pressurized 
(gEOS = 1-0), gas-rich major merger models with (model 
GE, right panel) and without (model GEn, left panel) su- 
permassive BH feedback. In both cases, stellar disks form 
from gaseous remnant disks created during the merger. 
The presence of supermassive BHs slightly decreases the 
mass of the central spheroid and increas es its effective ra¬ 
dius, i n agreement with the findings of IRobertson et id1 
il2nn^ for mergers that produce elliptical galaxies. Gom- 
bined, the velocity fields shown in Figure 0 and the stel- 
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Fig. 5. — Rotation curves of remnants from mergers of galaxies with moderately pressurized ISM and growing supermassive black 
holes for a variety of disk orientations. Shown is the rotation for new stars formed after the merger (open circles), old stars formed before 
the merger (open triangles), and all stars (solid circles). The formation of rapidly rotating remnants is not limited to prograde-prograde 
coplanar mergers (model GC, lower right panel), but can occur in both polar (model GCoF, upper right panel) and inclined (model GGoO, 
lower left panel) orbits. Gas-rich mergers can also produce remnants with very little rotation (model GCoE, upper middle panel) or 
counter-rotating cores (model GGoP, lower middle panel). 


lar disks shown in Figure^ldemonstrate that disk galaxy 
remnants containing rapidly rotating, stellar disk com¬ 
ponents can form from high-angular momentum, gas-rich 
mergers when the ISM is pressurized. 

3.2. Orbital Study 

While the results of 30 demonstrated that the pres¬ 
surization of the multiphase ISM from star formation 
and feedback can produce remnant disks in gas-rich, 
prograde-prograde coplanar mergers, the generality of a 
merger-driven scenario for disk galaxy formation would 
require a wider range of encounters to produce r emnant 
disks. Using simulations without star formation. iBarn^ 
ll2002^ demonstrated that gaseous disks could form in 
mergers that were not exactly coplanar, including polar 
orbits. Here, we examine the impact of varying a range of 
progenitor disk orientations and pericentric passage dis¬ 
tances on remnant disk formation including the effects of 
star formation and energetic feedback mechanisms. 

Figure El shows the rotation curves for six dif¬ 
ferent fairly radial encounters, including a prograde- 
prograde coplanar merger (model GC, lower right panel), 
retrograde-retrograde coplanar (model GCoC, upper left 
panel), polar (models GCoF, upper right panel, and 
GCoP, middle lower panel), and intermediate cases 
(models GCoE, upper middle panel, and GCoO, lower 
left panel). In each case, the same moderately- 
pressurized ISM ((/EOS = 0.5) is used and enables the 


new stellar component to rapidly rotate. Remarkably, 
some highly non-coplanar encounters produce remnants 
whose entire stellar remnant is rapidly rotating (models 
GCoF and GCoO) and demonstrate that merger-driven 
disk formation is not limited to the s i ngle c oplanar or¬ 
bit discussed by [SEnSge^^^ernguiil (|200^. Certainly, 
some orbits lead to remnants with almost no average ro¬ 
tation (model GCoE), though with a rapidly rotating 
young stellar component. Interestingly, the retrograde- 
retrograde polar encounter (model GCoP) produces a 
remnant with a new stellar component that counter¬ 
rotates in its interior relative to the outer young stel¬ 
lar material. Previous studies have demonstrated the 
formation of counter-rotating disks to form in rem- 
nants of mergers involving progenitors with /gas ss 0.1 
lIHernauist fc Bar nesll99 1i) or in spiral galaxie s via other 
mechanisms (e.sr. lThakar fc R.vdeiJfl99011199'^ . We plan 
to examine the creation of counter-rotating cores through 
extremely gas-rich major mergers in subsequent work, 
but we note here that the relative ages of the older, 
central disk in model GCoP compared with its outer, 
younger disk is consistent with the picture for the for- 
mation of kinematically d ecoupled cores presented by 
iHernauist fc Barn^ lll991i) where the central disk forms 
primarily from strongly-torqued material during the first 
passage and the outer disk forms later from material that 
retains some of its original orbital or disk angular mo¬ 
mentum. 



























10 


Robertson et al. 



-10 -5 0 5 -10 -5 0 5 10 -5 0 5 10 

X [h'* kpc] X kpc] x [h'* kpc] 


Fig. 6. — Rotational support V-cot/o' and projected stellar mass for three orbits that produce remnant disks in very gas-rich mergers. 
Shown are the logarithmic surface density of the stars (upper panels) and the Frot/cr measured for all stars (solid circles), and stars formed 
before (open triangles) or after (open circles) the merger of galaxies with moderately pressurized ISM models (^eos = 0.5). The polar (model 
GCoF, left panels) and inclined (model GCoO, middle panels) orbits lead to remnants that are disk-dominated and rotationally-supported. 
The prograde-prograde coplanar orbit (model GC, right panels) also produces a rotationally-supported disk system. 


The rotational support of the remnants GCoF and 
GCoO formed from non-coplanar mergers, shown in Fig¬ 
ure El with the coplanar encounter GC, is surprising but 
bolsters the merger-driven scenario for disk formation. 
In each merger, the total stellar rotational support ex¬ 
ceeds Vrot/o" = 1 and the new stellar component peaks 
at Vrot/CT > 5. While the mass of the prograde-prograde 
coplanar merger remnant is still dominated by the cen¬ 
tral bulge (bulge-to-disk ratio B : D 1.2 : 1.0), 
the other models are dts/c-dominated with bulge-to-disk 
ratios of R : H « 1.0 : 2.1 for model GCoF and 
B ■. D Ki 1.0 : 1.6 for model GCoO including thick 
disk components. At the time of the merger these non- 
coplanar mergers have /gas ~ 0.6 — 0.7 of their bary- 
onic component in gas, compared with /gas ~ 0.5 for the 
coplanar merger. The rotationally-supported remnants 
from non-coplanar orbits experience less star formation 
duri ng their first passage tha n the coplanar orbit (see 
also iMihos fc Hernauistlll996tl . leaving more gas avail¬ 
able to form a remnant disk after the final coalescence. 
The inclined merger model GCoE (Figure El upper mid¬ 
dle panel) that produces a remnant with very little net 
rotation has a /gas ~ 0.2 immediately after the final co¬ 
alescence, whereas the disk remnants formed by models 
GCoF and GCoO have /gas ~ 0.5 and /gas « 0.4, re¬ 
spectively. Clearly, both the amount of gas present in 
the interacting system before the height of merger and 
the fraction of gas consumed during the final coalescence 
may influence the formation of remnant disks by improv¬ 


ing the rotational support of the remnant. 

The simulations that explore the effects of disk orienta¬ 
tion demonstrate that for a given amount of orbital angu¬ 
lar momentum in a merger the contribution of disk angu¬ 
lar momentum can affect the presence of a remnant disk, 
both by providing angular momentum to the remnant 
gaseous disk that forms the rapidly-rotating stellar com¬ 
ponent and by influencing the amount of gas consumed 
during the merger. However, the orbital configurations 
considered above only alter the disk orientation and not 
the orbital angular momentum of the encounter. Chang¬ 
ing the pericentric passage distance can influence both 
the angular momentum of the remnant and the amount 
of gas available to form a remnant disk both by chang¬ 
ing the timescale of the merger and by affecting the gas 
consumption during the interaction. 

Figure [3 shows the rotation curves of remnants with 
moderately pressurized ISM models (gsos = 0.5) formed 
in prograde-prograde coplanar gas-rich (/gas = 0.99) 
mergers with pericentric passage distances increasing 
from Tperi — 7?d, progenitor “ 2.Th kpc (model GCrA, 
upper left panel) to Tperi — 0.37?vir,progenitor ~ 48/l 
kpc (model GCrD, lower right panel). Also shown for 
comparison is a polar merger with Tperi ~ 2i?d, progenitor 
(upper right panel), the same pericentric passage dis¬ 
tance as the prograde-prograde coplanar model GC (up¬ 
per middle panel). The nearly head-on collision (model 
GCrA) produces a remnant with almost no net rotation; 
after the height of the merger, only a small amount of gas 
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Fig. 7. — Rotation curves of remnants from mergers of galaxies with moderately pressurized ISM and growing supermassive black holes 
for a variety of pericentric passage distances. Shown is the rotation for new stars formed after the merger (open circles), old stars formed 
before the merger (open triangles), and all stars (solid circles) for pericentric passage distances increasing from rpe^; = = 2.7h~^ kpc 

(model GCrA, upper left panel) to rperi = O.SRvir = 48h“^ kpc (model GCrD, lower right panel). The formation of rapidly rotating 
remnants is appears to require orbits that are not too radial, since the first passage is violent and too much gas is converted to stars before 
the merger. Orbits also cannot be too wide, since the merger timescale is extended and results in a mostly stellar final merger. Remnants 
from mergers with pericentric passage distances of a few progenitor disk scale lengths, such as the coplanar (model GG, upper middle 
panel) and polar orbits shown here (model GGoF, upper right panel), appear to be most efficient at producing remnant disks. 


remains (/gas = 0.13). The gas available for a remnant 
disk increases to /gas ~ 0.5 — 0.6 when the pericentric 
passage disk scale length is roughly twice the disk scale 
length (models GC and GGoF), and decreases towards 
larger Tperi with /gas ~ 0.2 — 0.4 (models GGrB-GGrD). 
These orbits have a substantially longer merger timescale 
than the more radial orbits, allowing for more quiescent 
star formation and a larger stellar component present in 
the progenitors at the height of the merger. The com¬ 
peting effects of increasing the available orbital angular 
momentum by increasing rperi and decreasing the gas 
available to form a remnant disk owing to extended qui¬ 
escent star formation then lead to a range of rperi from 
which remnant disks are likely to form. While additional 
simulations will be needed to better characterize the best 
orbits for remnant disk formation, from the simulations 
performed here we estimate that fairly radial orbits with 
pericentric passage distances a few times larger than the 
disk scale length are most ideal for merger-driven disk 
formation in gas-rich mergers. 

3.3. Gas Fraction Study 

The preceding feedback and orbital studies have 
demonstrated that gas-rich mergers with a pressurized 
ISM can produce remnant disks for a variety of orbital 
configurations. The results of those simulations suggest 


that the primary influence on the presence of remnant 
disks is the amount of gaseous material left after the fi¬ 
nal coalescence; increased feedback or favorable (though 
diverse) orbits reduce the SFR during the first pericentric 
passage and the height of the merger, leaving enough gas 
after the merger to form a rapidly-rotating stellar com¬ 
ponent. By altering directly the gas fraction of the pro¬ 
genitor systems, the role of gas fraction on the rotational 
support of merger remnants can be determined. 

Figure IHl shows the impact of increasing the gas frac¬ 
tion from /gas = 0.4 (left panel, model DC) to /gas = 0.8 
(right panel, model FC) in the progenitor disks. Shown 
are prograde-prograde coplanar mergers with a moder¬ 
ately pressurized ISM model {qeos = 0.5). While the 
new stars in each model produce a rotationally-supported 
structure, the average rotational support of the system is 
well correlated with the gas fraction as judged from the 
increasing K-ot/o’ as a function of progenitor /gas. The 
gas fraction of the merging galaxies before and after the 
final coalescence are nearly in proportion to the original 
gas progenitor gas fractions, with /gas = 0.2 before and 
/gas = 0.1 after the final merger in model DC and roughly 
twice that those values for the FC model with double the 
original gas content. However, the original /gas = 0.8 
progenitors in the model FC merger are not gas-rich 
enough to produce a remnant that has a rotationally- 
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Fig. 8.— Rotational support Vrot/o' and projected stellar mass for mergers with progenitor gas fractions ranging from /gas = 0.4 (model 
DC, left panel) to /gas = 0.8 (model FC, right panel), mergers. Shown are the logarithmic surface density of the stars (upper panels) and 
the Vrot/cr measured for all stars (solid circles), and stars formed before (open triangles) or after (open circles) the merger of galaxies with 
moderately pressurized ISM models (q'eos = 0.5). The rotational component of these dispersion dominated galaxies correlate with the gas 
fraction, and progenitor galaxies will likely need an initial gas fraction of /gas > 0.8 to produce remnant disks (see Figure EJ. 


supported total stellar component with this orbital con¬ 
figuration. These simulations suggest that given the typ¬ 
ical ^ Gyr timescales involved, major merger progenitors 
must have original gas fractions /gas ^ 0.8 in order to sat¬ 
isfy the /gas ^0.5 requirement during final coalescence 
for a disk-dominated remnant. 

3.4. Minor Merger Study 

The simulations of i13.1l3.3l have demonstrated that 
gas-rich mergers can lead to remnant disk systems with 
substantial rotational support. However, these simula¬ 
tions have modeled only equal mass encounters. While 
major mergers on average contribute significantly to the 
final mass of an average galaxy at the present day, minor 
mergers between un equal mass galaxies are cosmologi- 
cally more frequent lILacev fc ColellTQ^ . 


The survival of disks in minor merger or infall 
scenarios have been readily addressed both semi- 


analvticallv JToth & Ostrikeil119921 

IBenson et alJ 

2004^ 

and numerically llOuinn et a,l .111 9931 1 Walker et all 

|l99fi|: 

iVelazQuez & Whitelll999l 

Font et al.ll2001l:IBarnesl 

2002^. 

Recently. iBournaud et al. 

1 I 2 OO 5 II has simulated compar- 


atively gas-poor (/gas < 0.2) mergers between galaxies 
with mass ratios in the range 4 : 1 — 10 : 1 which 
produce remnants with disk-like stellar structures bu t 
elliptical-like kinematics (see afso lRourna.ud et a, 1.11200^ . 
Currently, the combined results of these analyses sug¬ 
gest that stellar minor mergers with disks will pro¬ 
duce remnant systems that are kinematically dispersion- 


dominated but the cosmological infall of satellites ex¬ 
pected in the ACDM cosmology is consistent with the 
observed distribution of thicknesses for disk galaxies. 

Of additiona l interest to a scenario for merger-driven 
disk formation. iBarne i(| 2003 ) demonstrated in hydrody- 
namical simulations without star formation that gas-rich 
minor mergers can produce small gaseous disk remnants. 
A cosmological scenario for merger-driven disk formation 
would likely involve such minor mergers given their abun¬ 
dance in the mass-accretion histories of galaxies. Figure 
13 shows the rotational support of remnants formed in 
prograde-prograde coplanar mergers of the Wir = 160 
km s“^ galaxy model considered in H3.113.31 with the 
Wir = 80 km s“^ galaxy described in E3 These minor 
mergers with mass ratio m 2 : mi Ri 8 : 1 use a moder¬ 
ately pressurized ISM (^eos = 0.5) and progenitor gas 
fractions increasing from /gas = 0.4 (model DCm, left 
panel) to /gas = 0.8 (model FCm, right panel). The rem¬ 
nants are disk-dominated and rotationally-supported, 
with Vrot/o" ~ 2 at 5h~^ kpc. The smaller system is 
still massive enough to foment central star formation by 
driving gas inwards in the larger galaxy, increasing the 
velocity dispersion by a factor of « 2 to cr « 80 km s“^ at 
the center each remnant. These results are qualitatively 
consistent with previous work simulating the response 
of galaxy disks to minor mergers a t lower gas fractions 
(fZs ~ O-l. e.g. iHernauistl 119891 iMihos fc Hernauisd 
Il994bt iHernauist fc Mihc^ I1995I1 . but in our substan¬ 
tially more gas-rich systems the additional supply of gas 
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Fig. 9. — Rotational support Viot/o" and projected stellar mass for three minor mergers (mass ratio m 2 ■ mi =8:1) for progenitor 
gas fractions ranging from /gas = 0.4 (model DCm, left panel) to /gas = 0.8 (model FCm, right panel). Shown are the logarithmic surface 
density of the stars (upper panels) and the Dot/o' measured for all stars (solid circles), and stars formed before (open triangles) or after 
(open circles) the merger of galaxies with moderately pressurized ISM models (gsos = 0.5). The minor mergers produce disk-dominated, 
rapidly-rotating remnants with Dot /o' > 2. The central velocity dispersion of the larger galaxy roughly double during the merger, but 
the rotational suppor t of the remnants are considerably larger than the Drot/o measured for comparable gas-poor minor mergers (e.g. 
IBournaud et al.ll2n05ll . 


allo'ws the larger galaxy to retain comparatively more 
mass-'weighted stellar rotation by the end of the simu¬ 
lation through continued star formation in the extended 
gaseous disk. The velocity dispersion in the outer re¬ 
gions of our simulated remnants increases from ct ~ 30 
km s“^ in the larger progenitor at 5h~^ kpc before the 
merger to cr « 50 — 60 km in the remnant disk, de¬ 
creasing 'with progenitor gas fraction. While the remnant 
properties are not strong functions of the gas fraction for 
our very gas-rich mergers, with each remnant display¬ 
ing similar rotational support and bulge-to-disk ratios 
(B : I? ss 1 : 2 — 1 : 3, see Table |21l, the Trot/o’ « 2 — 2.5 
throughout most of the remnant disks substantially ex¬ 
ceeds the k/.Qt/o’ « 1.6 for the comparati'v el'V' gas-poor 
mergers reported by IBournaud et al.l 120051) for mergers 
of mass ratio m 2 : toi = 8 : 1 (see their Figure 10). 
We therefore infer that increasing the gas fraction of the 
disk progenitors of minor mergers b e 'vond the maximum 
/gas = 0.2 used by IBournaud et al.1 11200,51) will improve 
the rotational support of remnant disk galaxies, but we 
cannot fully detail the strength of the effect with the 
simulations presented here. However, based upon the 
success of the gas-rich minor mergers in retaining signifi¬ 
cant rotational support for the remnant disks we suggest 
that gas-rich minor mergers may be an arbiter of merger- 
driven disk formation and should be studied further. 


3.5. Small Progenitor Study 

The scale-dependent physics of gas cooling and star 
formation may influence merger-driven disk formation 
for a mass-sequence of galaxies. To examine this pos¬ 
sibility, we repeat the gas fraction study presented in 
H3.31 substituting smaller, Dvir = 80 galaxies for the pro¬ 
genitor systems while maintaining the same orbit (i.e. 
the same physical Cperi and Dnit). Figure [TOI shows the 
ViQtja as a function of radius for mergers with mod¬ 
erately pressurized ISM models (^eos = 0-5) and pro¬ 
genitor gas fractions increasing from /gas = 0.4 (model 
DCs, left panel) to /gas = 0.8 (model FCs, right panel). 
The remnants produced in these lower-mass mergers are 
rotationally-supported only in their outer regions. The 
merger with /gas = 0.4 progenitors (model DCs) does 
not contain a signihcant rotating component comprised 
from new stars, the only such system in our suite of sim¬ 
ulations. The higher gas fraction progenitors (models 
ECs and FCs) do produce rotating structures from their 
young stellar components, but the average rotation of 
these systems is dominated by rotation in their large 
spheroid components. 

As discovered for more massive progenitors in ii3.3l the 
gas fractions of the merging systems both before and af¬ 
ter the final coalescence scale roughly with the initial pro¬ 
genitor gas fraction for the same orbit and ISM pressur¬ 
ization. These systems each have gas fractions /gas < 0.5 
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Fig. 10. — Rotational support Vrot/o' and projected stellar mass for mergers with progenitor gas fractions ranging from /gas = 0.4 (model 
DC, left panel) to /gas = 0.8 (model FC, right panel), mergers for small V^ir = 80 km s~^ galaxies. Shown are the logarithmic surface 
density of the stars (upper panels) and the Vrot/o" measured for all stars (solid circles), and stars formed before (open triangles) or after 
(open circles) the merger of galaxies with moderately pressurized ISM models (^eos = 0.5). The rotational component of these dispersion 
dominated galaxies correlate with the gas fraction, though the relative importance of rotation is higher in these small mass systems than 
the more massive systems in similar mergers (see Figure . The progenitor galaxies will likely need an initial gas fraction of /gas > 0.8 
and the merging system will need /gas >0.5 just prior to the final coalescence to produce remnant disks (see Figure 13}. 


just before the final coalescence and /gas < 0.25 imme¬ 
diately afterwards, below the gas fractions measured in 
mergers that produce remnant disks. While the gas den¬ 
sities in the disks of these systems are lower than in the 
more massive systems considered in a and therefore 
produce relatively less star formation in quiescence, the 
merger timescale for the less massive systems is longer 
and leads to similar gas fractions at the final merger. If 
instead of using a physically identical orbit for both the 
lower-mass and higher-mass progenitors the orbits for the 
lower-mass progenitors were scaled with the disk scale 
lengths, or similarly the virial radii, the less massive sys¬ 
tems would have merged on a shorter timescale and may 
have produced more rotationally-supported remnants. 
Given the increased rotational support in the remnants 
of lower-mass progenitors relative to the higher-mass pro¬ 
genitors, we speculate here that lower-mass systems may 
be more conducive to merger-driven disk formation than 
can be inferred from the simulations we have performed. 
We leave explorations of these suspicions to future work. 

4. DISCUSSION 

The suite of mergers simulated in this work indicates 
that merger-driven disk formation may viable in the gas- 
rich mergers that build high-redshift galaxies. The re¬ 
quirements that the merging system is gas-dominated 
(/gas ^ 0.5) during the final coalescence to provide 


enough gas immediately after the merger to sustain a 
remnant disk is likely a necessary but not sufficient con¬ 
dition for merger-driven disk formation. However, phys¬ 
ical effects that increase the gas fraction during the final 
merger also improve the structure and rotational sup¬ 
port of remnant disks. Increased pressurization of the 
ISM, favorable orbits with less-violent first passages or 
shorter merger timescales, small mass ratio mergers that 
decrease the interaction-induced star formation in the 
more massive progenitors, and simply increased gas frac¬ 
tions in the progenitors all improve the rotational sup¬ 
port of the remnants and in some cases lead to disk- 
dominated, rapidly rotating stellar remnants. 

The importance of gas-rich mergers for a variety of 
galaxy properties h as been increasing ly realized (for a 
discussion see e.g., iBrook et all 1200,^ . and the results 
of our simulations of gas-rich progenitors bear on hi- 
erar chical model s for th e mass assembly of disk galax¬ 
ies. lAbadi et al.l ll2003bll report that a substantial frac¬ 
tion of the thick disk stars of a simulated disk galaxy 
formed in a cosmological setting originate as tidal de¬ 
bris fr om satellite systems accreted after the thin disk 
forms. IBrook et al.112004l2C)fiMl200(ill present a different 
scenario where the thick disk forms from multiple early 
gas-rich mergers (/gas ~ 0.5; C. Brook, private commu¬ 
nication) and the thin disk forms later from the gaseous 
debris of these accretion events. A recent photometric 
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survey of t he thick and thin components of edge-on disk 
galaxies hv lYoanhim Kr. Dalea.ntonl and kinematic 

decompositions of the thick and thin di s k com ponents 
of two galaxies bv lYoachim fc DalcantonI ll2005D suggest 
that thick disks may be kinematically disjoint, with lower 
rotational velocities than galactic thin disks. This ob¬ 
servational evidence s upports the t hick di sk accretion 
scenario forwarded bv lAbadi et al.l ll2003bH . However, 
if substantial thick disk mass accreted after the thin disk 
becomes predominately stel lar then heating o f the thin 
disk will occur anyway le.g. lQiiinn et al' IliM . Instead, 
if the epoch of thick disk accretion occurs when the large 
thin disk is gas-rich, then our simulations of minor merg¬ 
ers suggest that the thick disk may be safely accreted 
without unreasonably increasing the thin disk velocity 
dispersion. 

Our study of the formation of remnant disks in 
gas rich mergers provides a compliment to studies of 
elliptical galaxy form ation in gas-rich merging (e.g. 
IRobertson et al.1 l200fili . Major mergers between disk 
galaxies with progenitor gas fractions /gas < 0.8 al¬ 
most uniformly produce spheroid-dominated systems 
as the systems are not gas-dominated at the time 
of final coalescence. Such spheroid-dominated rem¬ 
nants are simil ar to our lowe r gas fraction mod¬ 
els EC-FC. As IRobertson et al.l (l2006|i demonstrate, 
these spheroid dominated systems obey elliptical galaxy 
scaling relations such as the Fundament al Plane 
ijDiorgovski fc PavislITo^ iDressler et'ZIllQSlt). 

Moreover, other apparently separate considerations 
make the case that if ellipticals did form through merg¬ 
ers, the progenitors must have been gas-rich (though not 
gas-dominated as the simulations here show that these 
events m ay leave behind system s with large disks). For 
example, iHernouist et alJ lll993ll showed that gas frac¬ 
tions larger than or of order 20—30% are required if merg¬ 
ers between disk galaxies are to account for the relatively 
high phase space densities of local ellipticals. Moreover, 
this same criterion ensures t hat the remnants obey ob¬ 
served scaling relations (e.g. USE e^ggi^^lJ]2QQ0ji and 
match kinematic properties (e.g. iCox et aT]l2flflfili of el¬ 
lipticals. 

The even more gas-rich mergers with progenitor gas 
fractions /gas > 0.8 with disk-dominated stellar compo¬ 
nents will not satisfy these elliptical galaxy scaling re¬ 
lations owing to their rotational support. These gas- 
rich merging scenarios for elliptical galaxy formation 
and remnant disk galaxy formation are fully consistent, 
as the remnant stellar disks form in systems with gas 
fractions large enough to sustain remnant gaseous disks 
whereas elliptical galaxies are generated in highly dissi¬ 
pative events where the vast majority of the gas is con¬ 
sumed and a large gaseous remnant disk cannot form 
owing to the nearly complete depletion of gas. In gas- 
rich merger events that lead to elliptical galaxy forma¬ 
tion, almost all of the stars form before or during the 
final coalescence of the galaxies. The production of rem¬ 
nant disks in mergers requires that progenitors are com¬ 
parably even more gas-rich to enable the formation of 
a remnant gaseous disk massive enough to eventually 
form a rapidly-rotating stellar disk that can dominate 
the stellar structure of the remnant. Naturally, this re¬ 
quires that the gas fraction before the final coalescence 
is /gas > 0.5, in agreement with the results of this study. 


Effects that reduce the consumption of gas before the 
final coalescence assuage the merger-driven production 
of disks, and while our work concentrates on feedback 
and orbital effects as well as the gas fraction of the pro¬ 
genitors we note that additional structural properties of 
the progenitors that reduce gas consumption during the 
merger (e.g. bulges) may also play a role in the for¬ 
mation of remnant disks. Lastly, the range of Hubble 
types observed for spiral galaxies spans a considerable 
range in both rotational support (Vrot/o") and bulge-to- 
disk ratio. The merger-driven scenario for disk galaxy 
formation here may account for systems with bulge com¬ 
ponents but likely cannot explain very late-type systems 
(e.g. Sd galaxies), which presumably must form almost 
entirely through dissipative processes that do not pro¬ 
duce hot stellar components. 

We emphasize that feedback associated with black 
hole growth, while essential f or accounting for the o b- 
served properties of quasars llHopkins et al.1 l2005 dlallbl) 
and the red colors of e llipticals fe.g. lSpring-el et all2005al: 
iHopkins et aLll^Ohdi . has little impact on the global 
structural properties of the rem nants studied her e . This 
is c onsistent w i th th e work of IRobertson et alJ l|2flfl(Tfl 
and iCox et ;d] ll200(itl who found that black hole feed¬ 
back likewise has negligible consequence for ellipticals 
on scales of order the effective radius. The signa¬ 
tures of black hole growth on the stellar distribution 
of merger remnants thus appears subtle and may be 
restricted to modifying the characteristics of the cen- 
tral starburst population for med during the merger (e.g. 
iMihos fc Hernaui^ll994al^ . 

5. SUMMARY 

To address the often destructive impact of the fre¬ 
quent merging in the ACDM cosmology on spiral galax¬ 
ies, we propose a new, merger-driven scenario for 
disk galaxy formation at high redshifts that supple¬ 
ments the_sta£dur^_pictoebased_on_dissi 2 ational col¬ 
lapse ilWhite fc Reeslll9'7^ iBlumenthal et al.lll984^ . In 
this scenario, gas-rich mergers can form rotationally- 
supported gaseous structures from residual angular mo¬ 
mentum after the final coalescence. This rapidly-rotating 
material will form stars, and if enough material is avail¬ 
able to form a substantial rotationally-supported com¬ 
ponent the remnant can resemble a disk galaxy both 
structurally and kinematically. We perform a suite of 
twenty-eight merger simulations to address the feasibility 
of such a scenario for disk galaxy formation and demon¬ 
strate that remnant disks can form from sufficiently gas 
rich mergers with a variety of orbits, mass ratios, and 
physical models for the ISM or the inclusion of growing 
supermassive black holes. We provide a detailed sum¬ 
mary of these results below. 

• We verify the results of b y iHarn ii iliiol and 
iSpringel fc Hernouistl ll2005ll in demonstrating that 
gas-rich mergers can produce remnant disks. 
Nearly every gas-rich merger produces a rapidly- 
rotating structure from stars formed in remnant 
gaseous disks after the merger. Furthermore, we 
demonstrate that certain gas-rich mergers can pro¬ 
duce remnants whose entire stellar structure is 
rotationally-supported. A necessary condition for 
the formation of these disk-dominated remnants is 
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inferred to be a gas fraction of /gas > 0.5 just prior 
to the final coalescence, which we have referred to 
as being “gas-dominated”. Merger-driven disk for¬ 
mation is shown to occur under a variety of gas- 
dominated mergers and is not limited to e.g. ide¬ 
alized prograde-prograde coplanar mergers. Favor¬ 
able equal mass encounters for merger-driven disk 
formation include certain polar and inclined orbits 
with limited star formation during the first peri¬ 
centric passage. A successful merger-driven disk 
formation scenario would requires extremely gas- 
rich progenitors with gas fractions /gas > 0.8, and 
is therefore likely limited to high redshift galaxy 
assembly. 

• The rotational support of remnants is shown to 
correlate with the pressuriza tion of the ISM. Using 
the m ultiphase ISM model of lSoringel fc HernauistI 
(1200,111 modified to allow for an adjustable pres- 
suriza tion in the form of an effective equation of 
state l|Soringel et al.ll^05bll . the merger simula¬ 
tions demonstrate that pressurized ISM models 
limit the amount of gas consumed before the final 
coalescence of the progenitors. The reduced star 
formation increases the gaseous material available 
to form a rapidly-rotating remnant disk and im¬ 
proves the rotational support of the remnant sys¬ 
tem. Weakly pressurized ISM models either pro¬ 
duce unstable progenitors owing to the high gas 
fractions or have efficient star formation. In both 
cases the merging systems experience a mostly stel¬ 
lar collision that simply kinematically heats the ex¬ 
isting stars without producing a newly-formed disk 
to contribute to the rotation of the remnant. 

• Gas-rich minor mergers can produce disk galaxies 
with substantial rotational support (Idot/cr > 2). 
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